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Abstract: We point out an inconsistency in perturbative QCD predictions previously 
used for dijet azimuthal decorrelations for azimuthal angles of A(/)dijet < 27r/3 between the 
two jets. We show how the inconsistency arises and how the calculations can be modified 
to provide more accurate results that exhibit a smaller scale dependence and give a better 
description of the data than the inconsistent results. We also explain how the quality of 
the predictions strongly depends on a perceivedly minor detail in the definition of the dijet 
phase space and give recommendations for future measurements. 
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1 Introduction 

Measurements of dijet azimuthal decorrelations in hadron-hadron collisions provide a unique 
testing ground for the predictions of perturbative quantum chromodynamics (pQCD). The 
dijet azimuthal decorrelation studies the production rates of dijet events as a function of the 
azimuthal angular separation between the two jets in an event that define the dijet system, 
” </>jet 2 |- The measured quantity, labeled P in this article and originally 
proposed by the D0 collaboration [1], is the dijet differential cross section, dudijet/dA^dijet, 
normalized by the inclusive dijet cross section, cidijet, integrated over A(/)dijet: 

P = — -(1,1) 

<^dijet dA^dijet 

The range of kinematically accessible values in A(/)dijet is indicated in figure 1 for 
processes with final states of different jet multiplicities. In 2 —)• 2 processes Ai^dijet has 
always the largest possible value of A(/)dijet = (figure la). If A^dijet is significantly below 
TT, then the quantity P is probing hard 2 —)• 3 and 2 —)• 4 processes, i.e. three-jet and four-jet 
production. Following the D0 measurement, the quantity P was also measured by the CMS 
and ATLAS collaborations [2, 3]. In all measurements the data are fairly well described 
by the theory predictions at next-to-leading order (NLO) pQCD for 37r/4 < A(^dijet < 

For smaller A^dijet, in particular for A(/)dijet < 27r/3, the theory predictions exhibit a large 
renormalization scale dependence and lie significantly below the data. 

In this article, we focus on the comparison of fixed-order pQCD predictions and data 
in the kinematic region of A()>dijet < 27r/3. In section 2 we introduce and compare the 
phase space definitions in the different analyses and discuss their effects on the kinematic 
constraints in 2 —)• 3 processes. In section 3 we show that the pQCD calculations by two 
of the experimental collaborations [1, 2] for the region of A(/>dijet < 27r/3 are inconsistent, 
and demonstrate how a correct treatment provides pQCD predictions with a reduced scale 
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Parameter 

Experiment (n 
D0 {pp, 1.96 TeV ) 

auction and center-c 
CMS {pp, 7 TeV) 

)f-mass energy) 
ATLAS {pp, 7 TeV) 

jet algorithm 
jet radius 

yinitial 

yfinal 

PTmin 

PTmax ranges 

Run II cone 

.^^cone = 0.7 

oo 

0.5 

40 GeV 

75-100 GeV 

100-130 GeV 

130-180 GeV 

>180 GeV 

anti-kj 

R = 0.5 

5.0 

1.1 

30 GeV 

80-110GeV 

110-140GeV 

140-200 GeV 

200-300 GeV 

>300 GeV 

anti-kt 

R = 0.6 

2.8 

0.8 

100 GeV 

110-160GeV 

160-210 GeV 

210-260 GeV 

260-310 GeV 

310-400 GeV 

400-500 GeV 

500-600 GeV 

600-800 GeV 

>800 GeV 


Table 1. Summary of the parameters defining the dijet phase space in the D0, CMS, and ATLAS 
measurements of dijet azimuthal decorrelations [1-3]. Variables are defined in the text. 


dependence. The results of these calculations also give a better description of the experi¬ 
mental data, as shown in section 4. In section 5 we discuss how a particular choice in the 
selection of the dijet phase space in the third experimental analysis [3] renders fixed-order 
pQCD predictions less accurate and how this can be improved in future measurements by 
a small modification in the dijet phase space definition. 


2 Phase space and kinematic constraints 

For a given process (e.g. pp or pp collisions) and center-of-mass energy, the measured quan¬ 
tity P, defined in equation (1.1), depends on additional choices, including the jet algorithm 
with its parameters, and the requirements on the jet rapidities y and the transverse jet 
momenta pT with respect to the beam direction. The initial jet selection may be carried 
out in a limited y region, with |y| < ^initial (where yinitiai can be adapted to the detector 
acceptance). The dijet system is then defined by the two jets with the highest px inside 
this region; here, these are labeled “jetl” and “jet2”. The final phase space for the ra¬ 
pidities yi ^2 of jetl and jet2 is then further constrained by |yi, 2 | < yfinai- Furthermore, the 
Pt of jet2 is required to be above a given threshold, PTmin, and the analysis results are 
presented in different regions of the px of jetl, pxmax- An overview of the choices for these 
parameters in the analyses by the D0, CMS, and ATLAS experiments is given in table 1. 
The main difference between the three scenarios regarding the scope of this article is the 
choice of yinitial- la the D0 scenario, the y region for the initial jet selection is unlimited 
(yinitiai = cc), while the ATLAS and CMS scenarios are limited to yinitiai = 2.8 and 5.0, 
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b) 2^3 



c) 2^4 
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Figure 1. Sketches of the azimuthal angular separation A(/)dijet between the two jets leading in px 
in an event for 2 —?► 2, 2 —>■ 3, and 2—^4 processes. Also indicated is the kinematically accessible 
range in Ac^dijet for the three configurations. 


respectively [3, 4], As a consequence of the choices for ^initial and PTmim the three scenarios 
then have different kinematic constraints for 2 —)• 3 processes as explained below: 

• Kinematic constraints for an nnlimited y region, yinitiai = oo 

For ^initial = CO, the selected jets, jetl and jet2, are always the two jets leading in px 
of the entire event. This selection criterion results in the kinematic constraint that 
the smallest possible A^dyet value in a 2 —>• 3 process (i.e. in a three-jet final state) 
is A^dijet = 27r/3 (cf. figure 1 b), while angles of Ai;f)dijet < 27r/3 are only accessible 
in final states with four or more jets (cf. figure 1 c).^ Therefore, for ^initial = oo, the 
dijet cross section for A(/)dijet < 27r/3 is a four-jet quantity, meaning that the lowest 
order pQCD contributions are from the four-jet tree-level matrix elements. 

• Kinematic constraints for a limited y region, ^initial < oc> 

If the y region for the initial jet selection is limited, it is possible that the two jets, 
selected for the dijet system, are not the two jets leading in px of the whole event. 
Table 2 gives an example for the ATLAS scenario, in which the leading jet in the 
event has \y\ > ^initial- In this case, the dijet system is made of the second and third 
leading jets, which are the two highest px jets inside the limited y region. Since there 
is no kinematic constraint for the azimuthal angular separation between the second 
and third leading jet, the region A(/)dijet < 27r/3 is also populated by three-jet final 
states. If such configurations are not prohibited by other phase space constraints, 
the dijet cross section for A(/)dijet < 27r/3 is a three-jet quantity. 

It depends on the requirements on yinitiai) (PTmax/\/s); and (pxmin/\/s)) whether a 
leading jet is kinematically allowed outside the region \y\ < ^initial and, as a consequence, 
three-jet configurations can populate the region of A(/)dijet < 27r/3. This can be tested 

^An event with exactly three jets can have Ai^dijet = 27r/3 only in a “Mercedes Star” configuration, 
where the jets have pxi = Pt2 ~ Pt3 and A(f)i^2 = At/ii^a = A(^2.3 = 27r/3. If the two jets leading in px in 
a three-jet event (with pxi > Px2 > Pxa) had A(()dijet < 27r/3, the vector sum of their transverse momenta 
could only be balanced, if the third jet had px3 > Px2, which would, however, contradict the assumption 
that px2 > Px 3 . 
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A 02,3 (radians) 
Ma’jet (TeV) 

0.528 • vr 

2.745 

xi (for Vs = 7 TeV) 
X 2 (for Vs = 7 TeV) 

0.990 

0.155 



jet 1 

jet 2 

jet 3 

Px (GeV) 

y 

4> (radians) 

405 

2.805 

0.000 • TT 

401 

-0.75 

0.920 • TT 

101 

-0.75 

1.448 • TT 


Table 2. The topology of an exclusive three-jet event, with the jet variables px, y, and (j> (left) and 
the event quantities A(^ 2 , 3 , three-jet invariant mass and the momentum fractions Xi and X 2 

for a center-of-mass energy of ^/s = 7TeV. In this event, the highest px jet is produced at large 
rapidity. If the dijet selection is restricted to jets with \y\ < j/initiai = 2.8 (as in the ATLAS scenario, 
see text), the selected dijet system does not include the highest px jet. This enables the azimuthal 
angular separation of the jets in the dijet system (here, A^dijet is determined by the azimuthal 
angle between the second and the third jet, A(() 2 , 3 ) to fall below the limit of Ac^dijet = 27r/3. 


using a cross section calculation based on tree-level 2 —>• 3 matrix elements as e.g. in NLO- 
Jet-| —h [5, 6]. We have used NLOJetH— h to compute the dijet differential cross section 
dcJdijet/dA^dijet for all three scenarios. The results for the ATLAS scenario are shown in fig¬ 
ure 2 and it is observed that up to and including the pxmax region of 400-500 GeV, the dijet 
differential cross section d(Tdijet/dA(?i)dijet receives non-zero contributions at A^dyet < 27r/3 
from three-jet final states. Therefore, in the ATLAS scenario, dudijet/dA^dijet is a three-jet 
quantity for all A(/)dijet in the pxmax regions with pxmax < 500 GeV. Only in the higher 
PTmax regions it becomes a four-jet quantity. In those regions, however, ATLAS has not 
published any measurement for A(/)dijet < 27r/3. 

Like ATLAS, the GMS scenario also has a limited y region for the initial jet selection, 
with lower requirements for pxmax and PTmiru but with a larger value of yinitiai = 5.0. We 
have computed dudijet/dA^dijet for the GMS scenario as well and find that in all pxmax 
regions the 2 —)• 3 tree-level predictions for dudijet/dAc/idijet are zero for A^dijet < 27r/3. In 
other words, in both the GMS and the D0 scenarios dudijet/dA^dijet is a four-jet quantity 
for A(/>dijet < 27r/3. 

We summarize our findings as follows: 

• The denominator of P, Udijet) is the inclusive dijet cross section, which is a two-jet 
quantity in all scenarios. 


• For A(^dijet > 27r/3, the numerator of P, diTdijet/dA^dijet) is a three-jet quantity in 
all scenarios. 

• For A(^dijet < 27r/3, the numerator of P is a four-jet quantity, if the initial y region 
is unlimited (yinitiai = 00 ) as in the D0 scenario, or if the yinitiai and px requirements 
prohibit the two jets with the highest px’s in an event from having \y\ > Pinitiah as 
in the CMS scenario. 

• If the Pinitiai and px requirements allow one of the two jets leading in px to have 
\y\ > ^initial) then the numerator of P is a three-jet quantity for all A(/)dijet- This is 
the case in the ATLAS scenario for the pxmax regions up to 400-500 GeV in pxmax- 
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Figure 2. The pQCD predictions of order 0{al) for the dijet differential cross section dcr/dAt/fdijet, 
as a function of in different regions of pTmax for all analysis bins of the ATLAS measurement. 

The hgure demonstrates that the 0{al) contributions to bins with A(/)dijet < 27r/3 and pxmax < 
500 GeV are small but non-zero. 

3 Perturbative QCD calculations for cross section ratios 

The pQCD prediction for a ratio R of two cross sections a a and <7^ in a given relative order 
of as (e.g. LO or NLO) can be computed from the ratio of the pQCD predictions for aA 
and as- For this purpose, both must be computed at the same relative order, which is not 
necessarily the same absolute order in A LO result is then given by iixo = 
and a NLO result by i?NLO = If numerator and denominator are calculated 

in different relative orders, cancellation effects between theoretical uncertainties can be 
compromised leading to an artificially increased renormalization scale dependence of the 
results as discussed with respect to jet shapes in sections 3.1 and 4 of reference [7]. 

For two-jet quantities, the LO and NLO pQCD predictions are given by calculations 
to order 0(0^) and 0(af), respectively. For each additional jet required for the final 
state, the respective powers of as increase by one, so that for example the LO (NLO) 
predictions for three-jet quantities are given by pQCD calculations to order 0{al) (0(af)). 
Combined with the findings from section 2, we obtain the rules for the calculation of the 
LO and NLO results for the quantity P in the three scenarios and in the different regions 
of Ai;i)dijet- These rules are listed in table 3 and compared to the computational procedures 
applied in the experimental publications [1-3]. The theory results published by DO and 
CMS for A(/)dijet < 27r/3 and labeled “NLO” in references [1, 2] are inconsistent, because 
they mix relative orders for the numerator (LO) and denominator (NLO). Replacing the 
NLO result for the denominator (in 0{a^)) by the corresponding LO (0(a^)) provides 
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scenario 

Aipdijet range 

order for P 

numerator 

denominator 

used in publication 

ATLAS 

all Af^dijet 

LO 0{as) 


0{a^s) 

the published LO and 

NLO are correct 

NLO 0{al) 


0{al) 

D0 

and 

CMS 

> 27r/3 

LO 0{as) 

0{al) 

0{a^s) 

the published LO and 

NLO are correct 

NLO 0{al) 

0{at) 

Oial) 

< 27r/3 

LO Oia^s) 

0(at) 

0{a^s) 

numerator: LO 0{a'^) 
denominator: NLO 0{al) 
(inconsistent, using 
mixed relative orders) 

NLO 0{al) 

Oial) 

Oial) 


Table 3. Correspondence between absolute orders in as in the calculations of numerator and de¬ 
nominator and the relative order in the quantity P. The right column comments on the calculations 
used in the experimental publications. 


the correct LO result for P below Ai;/)dijet = 27r/3. Alternatively, the correct NLO results 
at A(/)dijet < 27r/3 can be obtained by replacing the four-jet LO (0 (q!^)) results by results 
based on the four-jet matrix elements at NLO pQCD {0{a^)), which have become available 
in the last years [8, 9]. 

4 Results 

Following the prescriptions in table 3 we have computed the LO and NLO pQCD predic¬ 
tions for P in the DO, CMS, and ATLAS scenarios in the different A(/)cjijet regions. For 
comparison, we also derive the inconsistent “mixed-order” results for P as published by 
DO and CMS. 

All calculations are made in the MS-scheme [10] and for hve massless quark flavors, us¬ 
ing NLOJetH— h [5, 6] interfaced to fastNLO [11, 12]. The results are obtained for renor¬ 
malization and factorization scales of fiR = = PTmax; with the MSTW2008NLO [13] 

parameterization of the parton distribution functions of the proton, and with ag evolved 
from a value of as{Mz) = 0.120 according to the two-loop solution of the renormaliza¬ 
tion group equation. The uncertainty due to the scale dependence is computed from the 
variations of the ratio P for correlated variations of the scales in the numerator and de¬ 
nominator of = hf = PTmax/2 and Hr = hf = 2pTmax- The ATLAS collaboration 
has published non-perturbative corrections [14, 15], which are applied to the pQCD results 
to get the final theory prediction. These corrections are typically below 1% and never 
larger than 3%. The D0 and CMS collaborations have not provided non-perturbative 
corrections. In these cases, the pQCD results are directly compared to the data.^ The the¬ 
oretical calculations in this study differ slightly from the calculations used in the CMS and 
D0 publications due to different choices of the parton distribution functions and as{Mz)- 
Furthermore, the D0 collaboration chose different renormalization and factorization scales 

^In reference [16] non-perturbative corrections for the D0 results are shown to be typically below 2% 
and never larger than 4%. In the CMS publication [2] the non-perturbative corrections are quoted to vary 
between —13% at Ai^dijet = tt/S and +4% at Ac^jijet = rr. 
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Dijet Azimuthal Decorrelations at Hadron Colliders 
D0 Vs = 1.96 TeV CMS Vs = 7 TeV ATLAS Vs = 7 TeV 



Figure 3. Measurements of dijet azimuthal decorrelations at hadron colliders from the D0, CMS, 
and ATLAS experiments (from left to right) are displayed as a function of the azimuthal opening 
angle A^idijet of the dijet system for different requirements of the leading jet px (different markers). 
The measurements are compared to theoretical predictions based on NLO (solid lines) or LO pQCD 
(dashed lines), depending on whether the measured quantity is a three-jet or four-jet variable, 
respectively. The scale dependence of the pQCD calculation is indicated by the shaded areas. 


of fiR = fip = PTmax/2, and the CMS collaboration applied non-perturbative corrections. 
For the purpose of the following discussion, these differences are negligible. 

The experimental results from the D0, CMS, and ATLAS measurements are displayed 
in figure 3 over the entire A^^ijet range. The data are compared to theory at NLO or LO, 
depending on the A(/)ciijet range and the scenario. Over the whole range of PTmax and 
A</>dijet) the theoretical predictions are in agreement with the data, except for the ATLAS 
data at small Ai;^dijet- 
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Figure 4. Ratios of data from different experiments (columns) to fixed-order predictions as a 
function of Ac^dijet, from low pTmax (bottom) to high pTmax (top). The ratios are shown in different 
regions of A^dijet for the pQCD predictions at NLO (open circles) and LO (full circles), and also 
for the case of mixing different orders in numerator and denominator (triangles) for A(/)dijet < 27r/3. 
For better visibility the full circles have been slightly shifted towards smaller values of A())dijet. The 
scale dependence of the different pQCD calculations is indicated by the corresponding lines. 


The region of small A(/)dijet) including the transition at A(/)dijet = 27r/3 and the effects of 
the inconsistent mixed-order predictions, are further investigated in the following. Figure 4 
shows the ratios of data over the different theory predictions for A(/)dijet ^ Svr/d. The 
ratios are computed for the NLO results, the LO results, and the inconsistent results from 
mixed relative orders. Also shown are the uncertainty bands due to the scale dependence 
of the different theoretical calculations. For A(/)dijet > 27r/3, in all scenarios the NLO 
pQCD predictions are compared to the data. For A(/)dijet > 37r/4, these give a good 
description of the data within scale uncertainties, which are below 5-10%. In the range 
27r/3 < Ai;/)dijet < 37r/4, the O(a^) (i.e. three-jet NLO) calculation for the numerator is 
running out of phase space for three-jet final states as Ac/idyet —S' 2tt/3. This causes the 
0(af) calculation to effectively become a four-jet LO calculation. In this A(/>dijet range the 
NLO prediction still describes the data, but with an increasing scale dependence of up to 


- 8 - 
























































































30% as A(/)dijet ^ 27r/3. 

For the CMS and D0 scenarios at Aijidijet < 27r/3, we first focus on the inconsistent 
mixed-order calculations as published by the experiments. Figure 4 shows that over most 
of the range (and in particular towards lower Ac/idyet) these predictions are significantly 
below the data even outside their large scale dependence, and they do not describe the 
A(/)dijet dependence of the data. Compared to the inconsistent mixed-order calculations, 
the correct LO predictions have a significantly reduced scale dependence, and they give 
a much better description of the data. While they still do not reproduce the A^dijet 
dependence, almost all individual data points agree with the LO prediction within the 
reduced scale uncertainty. 

Although, for the ATLAS scenario the pQCD predictions for A(/>dijet < 27r/3 are tech¬ 
nically still of NLO, their scale dependence is as large as that of the mixed-order predictions 
for the CMS scenario, and the description of the data by both are equally poor. 

5 Recommendations for future measurements 

In section 3 we pointed out that for the numerator of P in the ATLAS scenario the three- 
jet NLO cross section calculations formally are of NLO also for A(/>dijet < 27r/3. The 
results presented in section 4, however, demonstrate that these NLO predictions exhibit a 
larger scale dependence and that they give a worse description of the data than the LO 
predictions for the DO and CMS results. The difference between the ATLAS and the DO 
and CMS scenarios was traced back to the choice of yinitiai in the dijet selection as explained 
in section 2. In contrast to the DO and CMS scenarios, the kinematic constraints in the 
ATLAS scenario do allow 2 —)• 3 processes to give small, but non-zero contributions to the 
dijet cross section for A(^dijet < 27r/3. Therefore, in this A(/)dijet range, while formally being 
a NLO pQCD prediction, the 0{ag) calculation for the numerator effectively is only a LO 
prediction, since the 0{ag) terms contribute less than one percent. This “formally NLO but 
effectively LO” calculation for the numerator exhibits the typical large scale dependence 
of a LO calculation while the NLO predictions for the denominator have a reduced scale 
dependence, as typical for NLO calculations. As a consequence, the NLO prediction for 
the ratio P has a scale dependence, which is similar to that of the mixed-order calculations 
and larger than that of the LO predictions for the DO and CMS scenarios. 

Therefore, we strongly recommend that future measurements of dijet azimuthal decor¬ 
relations use values of ^initial that, together with the PTmin and PTmax requirements, do not 
leave any phase space for 2 —)• 3 processes below Ai^dijet = 27r/3. Technically, this can be 
investigated by using a phase space generator or a three-jet pQCD LO calculation for the 
numerator of P. 

6 Summary and conclusion 

Measurements of dijet azimuthal decorrelations at hadron colliders continue to be a testing 
ground for pQCD predictions at higher orders, beyond what is probed in inclusive jet 
and inclusive dijet production. In particular in the phase space region of Ac/idyet < 27r/3, 
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dijet azimuthal correlations are sensitive to the dynamics of final states with four or more 
jets. In all previous publications of azimuthal decorrelations, based on the quantity P = 
(l/<7dijet) • (dudijet/dA^dijet); this region was poorly described by theoretical predictions. 
In this article we have identified two reasons for this shortcoming. 

In the publications by D0 [1] and CMS [2] , the poor theoretical description of the data 
is related to the inconsistent mixing of different relative orders in as in the predictions for 
the ratio P. We have performed a consistent LO calculation by computing both, numerator 
and denominator, at LO. This correct LO pQCD prediction not only exhibits a smaller scale 
dependence, but also gives a better description of the experimental data for A(/)dijet < 27r/3. 

The improvement due to the consistent LO calculation can, however, only be achieved 
for dehnitions of the dijet phase space that ensure the two jets of the dijet system to be also 
the two leading px jets in the events. We strongly recommend for future measurements 
of dijet azimuthal decorrelations at small A(/>dijet to perform the initial dijet selection 
accordingly. 

If this is taken into account, the future usage of four-jet NLO calculations will provide 
NLO pQCD predictions for the whole A(^dijet range, extending precision phenomenology 
for dijet azimuthal decorrelations to the region A(/>dijet < 27r/3. Since in this A(^dijet 
region the quantity P is proportional to future measurements with higher statistical 
precision can also be used for novel as determinations. This recommendation also applies 
to measurements of dijet azimuthal decorrelations based on the quantity [17, 18] when 
this is measured for < 27r/3. 
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